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Abstract 

The universal properties and structure of halo nuclei composed of two neutrons (2n) and a 
core are investigated within an effective quantum mechanics framework. We construct an effective 
interaction potential that exploits the separation of scales in halo nuclei and treat the nucleus as 
an effective three-body system. The uncertainty from higher orders in the expansion is quantified 
through theoretical error bands. First, we investigate the possibility to observe excited Efimov 
states in 2n halo nuclei. Based on the experimental data, ^'^C is the only halo nucleus candidate 
to possibly have an Efimov excited state, with an energy less than 7 keV below the scattering 
threshold. Second, we study the structure of ^'^C and other 2n halo nuclei. In particular, we 
calculate their matter form factors, radii, and two-neutron opening angles. 
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I. INTRODUCTION 



There has been a considerable interest in physical systems with large scattering lengths 
recently. The scattering of two particles with short-range interactions at sufficiently low 
energy is determined by their S-wave scattering length a. If a is much larger than the range 
of the interaction tq, the system shows universal properties The simplest example is the 
existence of a shallow two-body bound state if a is large and positive, but there are many 
more, including the effects of a limit cycle [2| and the Efimov effect Q in the three-body 
system. 

The best known example of a nuclear system with a large scattering length is the two- 
nucleon (NN) system. There are two independent S-wave scattering lengths that govern 
the low-energy scattering of nucleons. Both scattering lengths are significantly larger than 
the range of the interaction tq ~ l/'ni,r ~ 1.4 fm, while the effective ranges are of the same 
order as tq. As a consequence, the description of few-nucleon systems in an expansion in 
ro/|a| is useful. It has successfully been applied to various two-, three-, and four-nucleon 
observables (See Refs. 0, S, [6| and references therein). 

Another type of system that can be described with these techniques are halo nuclei: a 
special class of nuclear systems which offer the possibility of exploring universal behavior 
0, B H, 10, m. Halo nuclei consist of a tightly bound core and a "halo" of lightly bound 



nucleons. They are characterized by their large nuclear radius compared to the radius of 
the core. Equivalently, the separation energy of the halo nucleons is small compared to 
the excitation energy of the core. This separation of scales allows for the use of effective 
theories, where one can assume the core to be a structureless particle, and treat the nucleus 
as a few-body system of the core and the valence nucleons. 

The ffist application of effective field theory methods to halo nuclei was carried out 



in Refs. |12l . [ISj, where the neutron-alpha system ("^He") was considered. More recent 
studies have focused on the consistent inclusion of the Coulomb interation in two-body halo 
nuclei such as the proton-alpha [l3] and alpha-alpha systems [l^. Three-body halo nuclei 
composed of a core and two valence neutrons are of particular interest due to the possibility 
of these systems to display the Efimov effect . Efimov found that in three-body systems of 
non-relativistic particles, if at least two of the three pairs of particles have a large scattering 
length \a\ compared to the range tq of the interaction, there occurs a sequence of three-body 
bound states whose binding energies are spaced geometrically between ^^/mrg and h'^ /ma?'. 
The number of bound states grows to infinity, with an accumulation point at the three-body 
scattering threshold, in the limit ±a — > oo. The sequence of three-body bound states have 
universal properties that are independent of the details of the two-body potential at short 
distances. The influence of long-range Coulomb interactions on the geometric bound state 
spectrum was recently investigated in a model study |16]. 

The ffist experimental evidence for Efimov states in ultracold Cs atoms has recently been 
found through their signature in three-body recombination rates [17]. This signature could 
be unravelled by varying the scattering length a over several orders of magnitude using a 
Feshbach resonance. For halo nuclei, the interaction strength can not easily be varied and 
one has to look for different signatures of the Efimov effect. Since the ground state of a halo 
nucleus can not be uniquely identified as an Efimov state for fixed a, it is customary to look 
for excited states that have the Efimov character. One can then consider a halo nucleus to 
display the Efimov effect if it has at least one excited state with universal properties. 

In this paper, we explore the occurance of the Efimov effect and its well known universal 
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properties for 2n halo nuclei with a core of spin zero. From the earliest studies of halo nuclei, 
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C has been suggested as a good candidate for Efimov states [18|, with future theoretical 



work supporting this prediction [19|, |20]. We critically examine these earlier studies and 
perform an improved analysis in the framework of an effective theory. The uncertainties of 
our leading order analysis are quantified through error bands. In the second part of this 
paper, we focus on the structure of 2n halo nuclei. In particular, we calculate their matter 
form factors, radii, and two-neutron opening angles. Finally, we end with conclusions and 
an outlook. 



II. THEORETICAL FRAMEWORK 

For our study, we choose the effective quantum mechanics framework of Refs. [g], 
which is equivalent to using field-theoretic language for the problem at hand. The short- 
range interactions characteristic of halo nuclei are then described using an effective interac- 
tion potential. The low-energy behavior of the system can then be reproduced with a level of 
accuracy proportional to powers of the low-momentum scale Miow over the high-momentum 
scale Mhigh- The theory is valid up to a momentum, Mhigh, at which the errors become of 
order one. For example, the two-body interactions of halo nuclei can be characterized by 
their large scattering lengths, a ~ 1/Miow, and their range, Tq ~ 1/Mhigh- Such systems 
need to leading order one coupling parameter, Cq, for each pair interaction tuned to repro- 
duce the scattering length. The range of the interaction enters at next-to-leading order. For 
a > 0, there is a two-body bound state, and the binding energy can then be found through 
the universal formula: 

where is the reduced mass. The dots indicate corrections of the order y/2fiB2/Mhigh ~ ro/a. 

For the large separation of scales involved in halo nuclei, zero-range interactions can be 
used in constructing the effective interaction potential. This leads to a separable potential 
made up of contact interactions in a momentum expansion. The 2-body S-wave potential 
to leading order is 

{p\Veff\ff)=Cogip)gip') + .... (2) 

where the dots indicate higher order momentum dependent interactions which we will ne- 
glect. In a future study, the errors could be systematically reduced by including the effective 
range correction to Eq. ([2]). g{p) is the regulator function (sometimes called the form factor) 
of the theory. Of course, the low-energy observables must be independent of the regular- 
ization scheme, and one can choose the scheme most convenient for calculations. We use 
a momentum cutoff scheme, multiplying the coupling parameter with a Gaussian regulator 
function 

g{p) = exp (^-j^^ . (3) 

where A is the cutoff parameter. This regulator function obviously suppresses the contribu- 
tions of momenta p,p' ^ A, where the effective potential would break down and no longer 
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be valid. A natural choice for the value of A is therefore A ~ Mhigh, but observables are 
independent of A after renormalization. 

This potential is then used in the solution of the three-body Faddeev equations in terms 
of the spectator functions Fi{q), which represent the dynamics of the core {i = c) and the 
halo neutron {i = n). To find the bound state of a halo nucleus composed of two valence 
neutrons and a core with spin zero, the resulting coupled integral equations are simply a 
generalization of the three-boson equation^ (see |2l|] and references within). We use units 
such that h = c = 1 and the nucleon mass m = 1: 
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where the shifted momenta tt, tti, and tt2 are given by: 
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In the above equations, i?3 > is the three-body bound state energy and A is the number 
of nucleons in the core. The free three-body propagators for a spectator neutron Gq and a 
spectator core Gq are: 

v4 + l 2 A + 2 
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(10) 



The effects of the interactions are contained in the two-body T-matrices which are obtained 
by solving the Lippmann-Schwinger equation for the neutron-neutron and neutron-core in- 
teraction with the effective potential, Eq. ([2]). In the kinematics of Eqs. (jUlS]), we have: 



tn{q';B3 



lA+1 

71 A 



exp (2/(aLA')) erfc ( v^/(|a„,|A) 



+ VK(g'; B-s) exp 2K(g ; B:,)/A' erfc \^2Eniq'; B^)/k 



^ In fact, the equations are the same for any bound three-body system of two identical particles and a core 
with spin zero, which interact through the pair- wise zero-range potentials given in Eq. 
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with 



and 
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-— exp (2/(aLA')) erfc (V2/(|a„„|A) 



+ V^c(g'; 53) exp 2E,(g'; E3)/A' erfc \^2E,{q'; B^)/K 



with 



EM':E)=B,^^q'\ 



(13) 



(14) 



where the index n, c indicates the spectator particle. In Eqs. (fTTl [T3l) . erfc(x) = 1 — 
(2/-\/7r) exp(— t^)rft denotes the complementary error function, which will go quickly to 1 
for X <^ 1. If the cutoff is chosen large compared to all momentum scales involved in the 

problem: A ^ \/e, these T-matrices reproduce the usual effective range expansion 

at leading order. 

A two-body bound state appears as a simple pole in the two-body t-matrix at an energy 
E = —B2. In our renormalization of the coupling parameter, we have tuned Co(A) so that it 
reproduces the scattering length needed to produce this two-body binding energy to leading 
order, as given in Eq. ([1]): 
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This renormalization works equally well for the case of virtual states, a < 0. Also, for 
A ^ the exp and erfc functions both quickly approach 1, and the resulting relations 

are analogous to Ref. [if. 

For most halo nuclei, the S-wave scattering length is not as well known as the two-body 
bound (virtual) state energy. Therefore, we will generally use the two-body energies, E^n 
and Enc, as input parameters, from which we can calculate the scattering length to leading 
order through Eq. ([1]), l/a„i = ±^/2fIniEni, where the -|- refers to a bound state and the — to 
a virtual state, i = noi c, and is the corresponding reduced mass. The difference between 
the two-body energy and the S-wave scattering length is higher order in the expansion in 
Miow/Mhigh ~ ro/\a\. 



III. UNIVERSAL PROPERTIES 



A. The Efimov effect in 2n halo nulcei 



The three-body binding energies are given by the values of B^ for which the coupled 
integral equations, Eqs. (jH [5]), have a nontrivial solution. In principle, Eqs. (jH [5]), should 
also include a three-body force term which is required for proper renormalization. However, 
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due to the limit cycle behavior of this three-body force, it is always possible to choose a 
cutoff where the three-body force vanishes. As a consequence, we can simply drop the three- 
body force and work with a finite cutoff A as in Ref . |2l| . Tuning this cutoff to reproduce a 
given three-body observable, we can predict other low-energy observables by using the same 
cutoff [H, 23|. The cutoff is inversely related to the interaction radius (see 24| and references 



within), thus taking the cutoff to infinity is equivalent to taking the range of our potential 
to zero. It is at this limit that the Thomas collapse jlll will occur and the energy of the 
three-body ground state will diverge. 

The Thomas collapse is closely related to the Efimov effect in that the deepest three- 
body bound states of the Thomas collapse can be identified with the deepest Efimov states 
[26| . The sequence of three-body Efimov states can be found from our integral equations 
with sufficiently large scattering lengths by finding the spectrum of binding energies for a 
fixed cutoff. By increasing the cutoff new three-body bound states appear in the spectrum 
at critical values of A, which are geometrically separated. Accordingly, the energies of the 
deeper bound states increase in magnitude. The Thomas effect is seen by the divergence of 
the deepest bound state energy for A — > cxd. However, the states below the natural cutoff, 
l/(mrQ), are physically irrelevant. They are outside the range of validity of the effective 
theory and can be ignored. 

Conversely, the sequence of three-body Efimov states have universal properties that are 
insensitive to the details of the two-body potential at short distances, and hence independent 
of the cutoff. One such property is that for the resonant limit, a — > ±00, at which there 
are infinitely many arbitrarily-shallow three-body bound states, the ratio of the binding 
energies of neighboring bound states approaches a universal factor, Aq, as the threshold is 
approached: 

Bf/B'^^^^^ — ^ (Ao)^ as n ^ +00 with a = ±00. (16) 

This universal scaling factor Aq depends on the masses of the particles. In our case, the 
masses of the two neutrons are equal, rui = m2, and the core mass A = m^/rrii. The values 
of B^^ for n = 1, 2, 3 as a function of the core mass A are shown in Fig. [T]for a finite value 
of A = 10. Note that because we have taken a ±00, A is given in units of an arbitrary 
momentum scale k. All energies are then given in units of k^. One interesting feature is 
the appearance of a minimum in the binding energy. The absence of this minimum in the 
i?3°'* curve is due to the fact that the magnitude of B^^ is near the order of magnitude of A, 
and details of the regularization scheme become important. Also, we see that the binding 
energy quickly reaches an asymptotic value for very large A, and diverges for A — > 0. 
The dependence of the discrete scaling factor, Aq, on the core mass A is well known (see, 
e.g. the review [ij. Fig. 52). We have checked our code by reproducing these results. In 
particular, the discrete scaling factor is largest for all equal masses, where it has the same 
value Ao ~ 22.7 as for three identical bosons. In the very large core limit, A ^ 1, the 
discrete scaling factor approaches 15.7. In the vanishing core limit, A <^ 1, it approaches 1 
as all three-body binding energies diverge. 



B. Possibility of Efimov excited states in 2n halo nuclei 

Our main aim in this section is to assess which halo nuclei have the possibility of possessing 
an excited Efimov state. The ground state energy and the two-body energies can not be 
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FIG. 1: Spectrum of three-body bound states, when two of the particles have equal mass mi = 772-2, 
as a function of the mass ratio, A = ms/mi, in the resonant limit a — > ±00. The cutoff parameter 
was fixed at a value of A = 10.0 (the units of A and -63"^ are arbitrary, for details see text). 

predicted by our theory and are taken from experiment. In other words, we would generally 
like to know what the values of the two-body energies must be, or correspondingly how large 
the scattering lengths must be, in order to produce an excited Efimov state, knowing the 
ground state binding energy. 



To this end, we construct the parametric region defined by the ratios, y End ^3 versus 



Y Enn/ B]P . The boundary curves representing the existence of an excited Efimov state for 
various values of the core mass are shown in Fig. [2J An analogous study was carried out 
in Ref. [l^ (see below). All points which lie within the boundary curve have at least one 
excited Efimov state above the state with energy -83"'*, while points outside the curve have 
no excited states above this state. The curve itself is built up of the points for which the 
bI^^^^ binding energy is equal to the scattering threshold; i.e. B^^~^^^ = Eni for Enn or Enc 
bound, where E^i is the larger of Enn and Enc, and B^^^^ = for Enn and Enc virtual. The 
boundary curves in Fig. [2] were found with n = 1 in order to minimize the regulator effects. 
Due to the scaling symmetry of the sequence of three-body bound states, the nth state can 
always be interpreted as the ground state and the {n + l)th state as the first excited state. ^ 
Here it is of interest to note that these results differ from the results found by Amorim, 
et al. in an analogous study [19I] using a hard momentum cutoff rather than the Gaussian 
regulator scheme. We found that the results agree almost exactly for a core mass equal to 
the nucleon mass, A = 1, but differ significantly for all other values of the core mass. While 
the qualitative conclusions on the likelihood of Efimov states in 2n halo nuclei are the same, 
the quantitative results are different. In fact, in doing the numerical calculations with a 
hard cutoff ourselves, we found results that match those presented here. 



^ Because of the regulator effects, the curves found with n — Q are slightly different from the curves in 
Fig. [21 The curves are practically the same for larger values of n, as the numerical values of -63"'' are 
much smaller than the cutoff A for n > 0. 





These results represent the leading order calculations with the effective potential de- 
scribed in Eq. ([2]). The theoretical uncertainty in calculating the binding energy is of the 
order ?s T'o/a, where tq is the effective range of the potential, and a is the scattering length. 
As stated before, the cutoff parameter is related to the inverse of the range of the potential, 
such that we can approximate tq ^ 1/A. However, it is important in this error estimate 
that the A used comes from the result with n = 0, the true ground state, rather than n = 1. 
This corresponds to taking the "natural" value for the cutoff A. We therefore estimate the 
leading fractional error of our boundary curves as ~ a/S/x^^T^^/A, for z = tt, and c respec- 
tively. The resulting boundary curve including leading order error bands, using the case of 
core mass A = 18 is shown in Fig. [31 This graph is a good representation of the error bands 
for other core masses. The uncertainty of our leading order calculation becomes large for 



values of y Enn/ B^'^ + Enc/ near 1 and greater. At this point the low-energy observ- 
ables approach the order of magnitude of A, where the effective potential is no longer a good 
description of the three-body system. 

Now we discuss the implications of Figs. [2] and [3] for the existence of excited Efimov states 
in halo nuclei in more detail. The four quadrants of the parametric plane in these figures 
correspond to the four different types of three-body halo nuclei, determined by the different 
types of two-body subsystems. The upper-right quadrant corresponds to both the n-n and 
the n-c subsystems being bound, and is accordingly called All Bound. The lower-right 
quadrant is that in which the n-n subsystem is bound, but the n-c subsystem is unbound, 
and receives the name Tango [23|. Of course, because we are concerned with 2n halo nuclei, 
where the n particle is truely a neutron, these two quadrants are not of much interest in the 
present study. The upper-left quadrant corresponds to the unbound n-n subsystem with a 
bound n-c subsystem, for which we use the name Samba as recommended in 28||. The final 
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FIG. 2: Boundary curves in the yEnc/B. 



vs. 



Enn/Blp^ plane, where the binding energy of 



the excited Efimov state B. 



(n+l) 



is exactly at threshold. Negative values on the axes correspond 



to virtual two-body states. Boundary curves shown for various core masses A =1, 9, 18, and 100. 
Experimental data shown for ^'^C, ^^C, ^^Li, ^^Be, and ^^Be are taken from Ref. 30|. 
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FIG. 3: Boundary curve in the y Enc/B^^ vs. y Enn/B^'^ plane with leading order error bands 
Boundary curve shown for a core mass of ^ = 18 with the experimental data for ^''C from Ref. 30] 



quadrant corresponds to the three-body systems for which none of the two-body subsystems 
is bound. This system is referred to as a Borromean system. 

We can now use our plot to analyze the likelihood of the Efimov effect for the different 
types of three-body systems, with a focus on 2n halo nuclei. As one would expect, the 
Borromean systems offer the smallest chance of having an excited Efimov state, as the two- 
body energy would have to be very small, or accordingly the scattering length very large, 
to produce even one excited state. However, this can be achieved in ultra-cold atoms, as 
the presence of so called Feshbach resonances allows one to tune the two-body scattering 



length to a very large value [29j. Interestingly, the Samba systems have the largest region 
supporting the occurance of excited Efimov states. As long as the n-c scattering length is 
large enough, there can be a large variation in the n-n scattering length that would still 
allow for the Efimov effect. This agrees with the findings of Efimov [31], that as long as 2 
of the 3 two-body interactions have a large scattering length, the sequence of three-body 
binding energies can occur. 

Looking at possible halo nuclei candidates, we have plotted the positions of ^°C, ^^C, 
"Li, ^^Be, and ^^Be in Fig.|2l using the experimental values of the " Nuclear Data Evaluation 
Project" of TUNL [sO] for the n-c and three-body ground state energy data, and the standard 
value of the n-n scattering length, a„„ = (— 18.7±0.6) fm 3l| to calculate the n-n two-body 
energy according to Eq. ([1]). The only halo candidate that has any possibility of an excited 
Efimov state is ^''C, due to the large uncertainty in the n-^^C bound state energy. We will 
return to this nucleus shortly. The positions and relatively small experimental errors in the 
other halo nuclei data rule out the chance of finding excited Efimov states in these nuclei. 
Other halo nuclei candidates which exist have values of the two-body energies which are too 
large to even appear on our plot. 
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FIG. 4: Binding energy of the C excited Efimov state as a function of the n- C bound state 
energy (sohd hne) with leading order error bands. The dashed hne represents the scattering 
threshold which is given by = Enc- The inset shows the excited state energy relative to the 
scattering threshold. 



C. Efimov excited state for ^'^C 



The central value for the n-^^C bound state energy, Enc = (162 ± 112) keV 32j, lies 



almost exactly on the boundary region for A = 18 in Fig. [3l The large error in this value, 
however, dips well into the region where at least one excited Efimov state can occur. The 
error in the three-body ground state energy of ^''C is small compared to E^c- Thus, we can 
calculate the value of the excited state energies as a function of Enc, using the standard 
value for a„„, and fixing our cutoff to reproduce the experimental value of the ground state 
energy B^^^ = 3506.0 keV [s^, 3^. The result is plotted in Fig. |H where the solid line 
is the excited state energy, and the dashed line represents the scattering threshold. The 
inset graph shows the excited state energy relative to the scattering threshold. We find only 
one excited Efimov state, existing when Enc < 165 keV. For larger values of Enc, the 
system moves outside of the boundary curve, and the excited Efimov state is destroyed. The 
binding energy relative to the scattering threshold is always below 7 keV, a value very small 
in comparison with the ground state energy. Also, the error bands are large compared to 
the relative energy of the exicted state to the scattering threshold, with the lower error band 
almost always below the scattering threshold. 

We have again estimated this error using the theoretical uncertainty of our effective 
potential. In this first order calculation, the uncertainty in binding energies calculated using 
the two-body effective potential of Eq. ([2]) is ~ ro/a. Our effective potential, made up of 
contact interactions, will break down for momenta of the order of the pion mass scale. We 
therefore use the inverse of the pion mass m.^ = 140 MeV to estimate the effective range 
ro ~ l/m^r. The uncertainty in the binding energy of the excited state is then the quadratic 
sum of the uncertainties from the n-n and n-^^C interactions: ^J2^nnEnn/fn'i + '^/•J-ncEnc/fn'i- 
These uncertainties are of the same order of magnitude as those found assuming that the 
effective range is related to the inverse of the cutoff, tq ~ 1/A. 
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This result is in overall qualitative agreement with the previous studies of Amorim et 



al. [19j and Mazumdar et al. [20[, who also found the presence of a very weakly bound 



excited Efimov state in for sufficiently small values of Enc- However, both of these 
studies have a larger value for the excited state energy, with the Mazumdar group also 
finding a second excited state for Enc < 100 keV. Also, the disagreement with the results 
of the Amorim paper mentioned before casts doubts on the quantitative results of 1^. as 
a more recent study from the same group |33[| suggests better agreement with the results 
presented here. This newer analysis is one of a few recent studies extending the trajectory of 
this excited state into the scattering region to explore the possibility of finding a resonance 



in the n-C scattering sector [33|, 13J, 135 



IV. FORM FACTORS AND MEAN SQUARE RADII 

We are now interested in calculating other low-energy physical properties of three-body 
halo nuclei, specifically the matter density form factors and the mean square radii. The 
information needed to calculate such quantities is held in the wave functions of the known 
bound states. In the previous sections, we only required the Faddeev spectator functions F„ 
and Fc. However, the full three-body wave function can be reconstructed from the Faddeev 
spectator functions (see Appendix |A] for details). Once the three-body wave function is 
known, the corresponding one- and two-body matter density form factors can be computed. 
Finally, the mean square distances between two of the three particles as well as the mean 
square distance of one of the particles from the center of mass can be extracted from the 
proper form factor. 



A. One- and two-body matter density form factors 

The three-body wave functions found in Appendix |X] can be used to calculate other 
low-energy properties of the three-body bound state. With the Jacobi momentum states 
it is straight forward to calculate the Fourier transform of the one- and two-body matter 
densities with respect to the momentum transfer squared. These are defined as the one- 
and two-body matter density form factors, Tiik"^) and Tniik"^) respectively, where i = n,c. 
In the wave functions, the p Jacobi momentum describes the relative momentum between 
the two particles in the chosen two-body subsystem, while q describes the momentum of the 
spectator particle relative to the center of mass of the two-body subsystem.^ Therefore, the 
one-body form factors can be obtained as follows: 

^,(A;2) = J d^p d?q ^,(p, q)m,{p, q- fc), (17) 

where i = n,c depending on the desired two-body subsystem. The two-body form factors 
can be solved similarly: 

^nc(A;') = / d'p £q ^n{p, g)*n(p - k, q), (18) 



^ Recall that we use the spectator notation for the wave functions, where the index i on refers to the 
spectator particle. 
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and 



^nn(fc') = J d^P d'q g)^c(p - k, q). (19) 

These relations can be simplified further by using the fact that at leading order only S- 
waves contribute. Consequently, we project the three-body wave functions onto the S-wave, 
and then perform the angular integrations. In our normalization, the wave functions then 
obey the relation: 

=47r< >, (20) 

where < ... > denotes the angular average. This relation is then substituted into the above 
form factor relations, and the trivial angular integrations can be performed. Furthermore, 
the form factors will be normalized in the end such that = 0) = 1, so any constant 

overall factor can be dropped. For the one-body form factors we have: 

J^i{k^) = j dp j dq q"^ j dx ^^(p, g)^i(p, y^q^ + - 2qkx), (21) 

and for the two-body form factors we have: 

^nc(A;') = Jdpp^ Jdqq' j dx ^nip, q)^niVP^ + k^-^pkx, q), (22) 



and 



y^nn{k'^) = j dpp^ j dqq^ j dx ^^(p, q)^cUp^ + P - 2pkx, q). 



(23) 



The expressions relating g) and \E'c(p, g) to the solutions -F„(g) and Fc{q) of Eqs. (jH 

[5]) are given in Appendix |X1 We are now in the position to calculate these form factors for 
halo nuclei. 

As a general example, we have plotted the form factors for the ground state of ^"C, in the 
low-energy region in Fig. [5l using a bound state energy of 3506.0 keV, a n-n two-body virtual 
energy of 116.04 keV, and a n-c bound state energy of 161.0 keV. The theory breaks down 
for momentum transfers of the order of the pion-mass squared (fc^ ~ 0.5 fm~^) where the 
one-pion exchange interaction cannot be approximated by short-range contact interactions 
anymore. Here one can see that for low momentum transfer the one-body neutron, jF„(fc^), 
and the two-body core-neutron, Tncik'^)-, form factors lie nearly on top of each other. This 
is due to the fact that the core consists of 18 nucleons and, therefore, the center of mass 
is very near the core. This fact is also seen in the shallow slope of the one-body core form 
factor, jFc(/c^). 

The theoretical error bands in the form factors are an estimate arising from the theoretical 
uncertainty of our two-body effective potential, Eq. ([2]). In this first order calculation, the 
uncertainty in the effective potential comes from the next term in the expansion, which is 
related to the effective range. Therefore, the theoretical uncertainty is ~ rg/a. As discussed 
above, we use the inverse of the pion mass m-Tr = 140 MeV to approximate the effective range, 
ro ~ l/rriT,. With the form factors normalized to = 0) = 1, and because \Enc\ > |-E„„|, 
the theory error is estimated to be ~ (1 — J-')y/2fj,ncEnc/^n- 
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0.7 








0.3 



FIG. 5: The various one- and two-body matter density form factors with leading order error bands 
for the ground state of ^'^C in the low-energy region: Tnn{k'^) black solid line; Tnc{k'^) black dotted 
line; J-n{k'^) lighter (red) dashed line; J-c{k'^) lighter (red) dot-dashed line. 

B. Mean square radii and geometry of 2n halo nuclei 

The mean square radii for our three-body bound states are calculated from the matter 
density form factors in the low momentum transfer region. The matter density form factor 
is defined as the Fourier transform of the matter density: 



with the normalization J-'{k'^ = 0) = 1. In the low momentum transfer region, the exponen- 
tial can be expanded, and assuming a spherically symmetric matter density, we see that the 
slope of the form factor determines the mean square radius (r^): 



Of course, the mean square radius acquired depends on the choice of one- or two-body 
form factor. Since p describes the relative momentum of the two particles in the two-body 
subsystem chosen, the slope of Tniik"^) will give the mean square distance between the two 
particles in the chosen two-body subsystem, either (r,^„) or (r^^). Analogously, because q 
describes the momentum of the spectator particle relative to the center of mass of the two- 
body subsystem, the slope of jFj(A;^), will give the mean square distance of the spectator 
particle from the center of mass of the two-body subsystem, either (r^_„„) or (r^_„c)- 
However, it is more useful to calculate the distance of the individual particles from the 
center of mass of the three-body bound state. If hi is the slope of the one-body form factor 




(24) 





(25) 
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FIG. 6: (a) The various radii of the three-body system, (b) Further geometry of the 2n halo 
nucleus, specifically looking at the two neutron opening angle 9nn- 



J-'i{k'^) at the limit k"^ = 0, the mean square radius of one of the bodies i 
three-body center of mass is given by: 



-66. 1 - 



rrii 



2m„ -|- rric 



n, c) from the 



(26) 



where rrii is the mass of the desired particle i, and m„ and rric are the neutron and core 
masses, respectively. The various radii of the three-body system are illustrated in Fig. [6](a). 

We have extracted the radii by fitting a polynomial in k"^ to the form factor results for 
small k"^. We have used polynomials of varying degree up to 5th order in k'^ in order to 
verify the stability and convergence of the fit. We have found a satisfactory stability in the 
slope when fitting to a polynomial to the fourth order in k'^, up to a value of fc^ at which 
the form factor has dropped less than 10 percent. 

As with the binding energies in the sequence of three-body Efimov states as discussed in 
Sec. nil At the mean square radii of these states also display universal properties that are 
insensitive to the details of the two-body potential at short distances. One such property is 
that for the resonant limit, a — > ±oo, at which there are infinitely many arbitrarily-shallow 
three-body bound states, the ratio of the radii of neighboring states approaches a universal 
factor as the threshold is approached. This universal scaling factor is exactly the inverse of 
the universal scaling factor found for the ratio of binding energies (see Eq. fll6p ): 



(n+l) 



(Ao 



1-2 



as n 



-oo with a = ±oo. 



(27) 



Therefore, we can construct a dimensionless quantity from the root of the product of the 
mean square radius and the three-body binding energy, a/ (r^) B^. The ratio of this quantity 
for neighboring states approaches unity in the resonant limit as the threshold is approached. 

The extracted radii for known halo nuclei are shown in Table [H As input we have used 
the standard value of the n-n scattering length, a„„ = (—18.7 ± 0.6) fm [3l[, to calculate 
the n-n two-body virtual energy, Enn = 116.04 keV, along with the experimental values 
of the n-c two-body energies, E^c, shown in the third column of Table [T] (negative values 
correspond to virtual energies). As a three-body input, the cutoff is tuned to reproduce 
the experimental ground state binding energy, Bf'\ shown in the second column of Table [B 
These experimental values for the two-body and three-body energies are taken from the 
most recent results of the "Nuclear Data Evaluation Project" of TUNL [sol], except where 
otherwise noted. In the last column, we have given the experimental values for the n-n 
mean square radius, as given by Marques et al. 3^. These experimental values for -\/(rf^ 
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were found using the three-body correlation study in the dissociation of two neutrons in 
halo nuclei, along with the two neutron correlation function. However, the large uncertainty 
in these values is indicative of the need for more precise measurements of the mean square 
distances in 2n halo nuclei. 

Our results agree overall with the study done by Yamashita et al. using a similar three- 
body model [s?!]. Our study expands on this previous work by showing the leading order 
theoretical uncertainty as well as the results for an excited Efimov state in the case of ^^C. 

The leading order theoretical error is again estimated by the uncertainty of the two-body 
effective potential, Eq. ([2]), which is ~ To/a, where tq is the effective range of the interaction, 
and a the scattering length. Using the inverse of the pion mass to estimate the effective 
range, tq ~ I/t/^tt, the uncertainty in the radii is then calculated from the greater of the 
error arising from the n-n or n-c interaction: \j2jjincEnc/m-K or \j2iJLnnEnn/rnn- 

We will now discuss the results for the various known three-body halo nuclei as shown in 
Table [H 

For ^^Li, there is a relatively large uncertainty in the experimental values of both the 
ground state energy, Bf^ = (247 ± 80) keV |30|, and the n-c virtual energy, with two 



competing values: E^c = (-25 ± 15) keV |30|, and E^c = (-800 ± 250) keV [38(. For 



this reason it is advantageous to plot the various mean square radii over the full range of 
potential Enc values. The results, using the central value for the three-body binding energy 
as input, B^^ = 247 keV, can be seen in Fig. [3, with error bands estimated from the 
theoretical uncertainty, as described above. These plots are a good general example of the 
relation between the mean square radii and the virtual two-body n-c energy for Borromean 
halo nuclei, where none of the two-body subsystems are bound. As the n-c virtual energy 



Nucleus 


^3 [keV] 


Enc [keV] 








VW) [fm] 




"Li 


247 


-25 


8.7±0.7 


7.1±0.5 


6.5±0.5 


l.OzbO.l 






247 


-800 [38] 


6.8±1.8 


5.9±1.5 


5.3±1.4 


0.9±0.2 


6.6±1.5 




320 


-800 [38] 


6.2±1.6 


5.3±1.4 


4.8±1.3 


0.8±0.2 






170 


-800 [38] 


7.9±2.1 


6.7±1.8 


6.0±1.6 


1.0±0.3 






1120 


-200 [39] 


4.1±0.5 


3.5±0.5 


3.2±0.4 


0.40±0.05 


5.4±1.0 


i2Be 


3673 


503 


3.0±0.6 


2.5±0.5 


2.3±0.5 


0.32±0.07 




18c 


4940 


731 


2.6±0.7 


2.2±0.6 


2.1±0.5 


0.18±0.05 




20c 


3506 


161 


2.8±0.3 


2.4±0.3 


2.3±0.3 


0.19±0.02 






3506 


60 


2.8±0.2 


2.3±0.2 


2.2±0.2 


0.18±0.01 






3506 


0.0 


2.7±0.2 


2.2±0.2 


2.1±0.2 


0.18±0.01 




20c* 


65.0±6.8 


60 


42±3 


38±3 


41±3 


2.2±0.2 




20c* 


1.02±0.08 


0.0 


130±10 


97±7 


93±7 


6.9±0.5 





TABLE I: Various mean square radii of different halo nuclei. The second two columns show the 
input values for the three-body ground state energy and the two-body n-c energy (negative values 
corresponding to virtual energies), respectively, as given by [sB], except where otherwise noted. 
The experimental values for the n-n root mean square radii, shown in the last column, are taken 
from [36]. The rows marked by ^Oc* show the results for the excited Efimov state of ^^C, with 
binding energy displayed in the second column, which is found above the ground state (-B3 = 3506 
keV). 
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FIG. 7: The various mean square radii for ^^Li as a function of the n-^Li two-body energy (negative 
values correspond to the virtual state) with error bands from the theoretical uncertainty. As input, 
the n-n two body energy Enn = —116.04 keV, and the three-body binding energy B^f' = 247 keV 
were used. 



decreases in magnitude the three-body bound state increases slowly in size, with a more rapid 
increase in size as the energy approaches zero and crosses over into the Samba configuration, 
where the n-c subsystem becomes bound. 

In Table [H we have highlighted, using the central value of the three-body binding energy, 
the central values of the competing n-c energies. While the two-body virtual energy reported 
in jssi, Enc = —800 keV leads to \/Jr^ = (6.8 ± 1.8) fm in close agreement with the 
experimental result ^y {fnn)exp ~ (^'^ -'- -'-•^) ^ definite conclusion can not be reached 
due to the large error bars of the radii. We have also listed the upper and lower limits of 
the three-body binding energy, 5f ^ = 170 and 320 keV, along with our preferred value of 
Enc = —800 keV, which shows a more general result: for halo nuclei, the larger the three- 
body binding energy, the smaller the mean square radii. In terms of the plots in Fig. [71 
using a larger (smaller) value for sf^ as input would shift the curve down (up) in each plot. 
Due to the large uncertainties in both the theoretical and experimental values for ^^Li, there 
exists a large range of Enc values which would produce a a/ {r^n) value in agreement with 



the experimental value of Marques et al. [36 



As another example of a Borromean halo nucleus, we see that the calculated result for the 
n-n mean square radius of ^^Be, \/lr^ = (4.1 ± 0.5) fm is smaller than the experimental 
value a/ {''^nn)exp ~ ^^"^ ^■'^) l'^^]^ t)ut still withiu one error bar. Here we have used 
the central value of the two-body n-c virtual energy as reported by [39], Enc = —200 keV. 
In using the resonant limit, Enc = 0.0 keV, we see that the largest theoretical value for 
a/ {rnn) = (4.6 ± 0.3) fm, which allows an unbound n-c two-body subsystem, is still smaller 
than the experimental value. Another reported value for the two-body virtual energy, Enc = 
(—1900 ± 500) keV [30|, would produce even smaller mean square radii, even farther away 
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FIG. 8: The various mean square radii for '^C as a function of the n-C two-body energy with 
error bands from the theoretical uncertainty. As input, the n-n two body energy S„„ = —116.04 
keV, and the three-body binding energy B^^ = 3506 keV were used. 



from the experimental value. 

We now turn our attention to the so called Samba halo nuclei in which the n-c subsystem 
is bound. As examples, we have listed the results for ^^Be and ^^C, using the central values 
of the experimental energies 30[ in Table [B 

Of greater interest are the results from the case of ^'^C, as the large uncertainty in the 



n-c energy, with two competing values, Enc = (162 ± 112) keV [30], and Enc = 530 keV [40 
suggests that we look at the mean square radii over a range of Enc values. The results, using 
the central value for the three-body binding energy as input, 5^°^ = 3506 keV, can be seen 
in Fig. [HI with error bands estimated from the theoretical uncertainty, as described above. 
These plots are a good general example of the relation between the mean square radii and 
the two-body n-c binding energy for Samba halo nuclei. As the n-c binding energy decreases 
in magnitude the three-body bound state decreases slowly in size, with a slightly more 
rapid decrease in size as the energy approaches zero and crosses over into the Borromean 
configuration, where the n-c subsystem becomes unbound. This suggests that as the two- 
body n-c state is more weakly bound, the particles must be closer together in order for the 
three-body state to be bound with the same energy. 

As was shown in Sec. IIII C\ there possibly exists one Efimov excited state in for 
Enc < 165 keV. The mean square radii for this excited state were calculated over a range of 
Enc values and plotted, with leading order error bands, in Fig. O Here we see the interesting 
phenomenon that although the radii remain relatively constant in the middle of the range of 
n-c energies which allow the excited state, as the endpoints are approached, the radii begin 
to increase rapidly, and then to diverge, both when Enc 165 and 0.0 keV. It is at these 
points that the system moves outside the boundary curve depicted in Fig. [3l the Efimov 
excited state is destroyed, and the three particles consequently fly apart. 
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FIG. 9: The various mean square radii for the Efimov excited state of ■^''C as a function of the 
n-^^C two-body energy with error bands from the theoretical uncertainty. As input, the n-n two 
body energy = —116.04 keV, and the three-body binding energy b!^^ = 3506 keV were used. 



In Table HI we have highlighted the result using the central value of the accepted n- 
c two-body binding energy, as well as two values which lead to an Efimov excited state, 
including the resonant limit E^c = 0.0 keV. The rows marked by ^°C* represent the results 
of the Efimov excited state. The three-body binding energy of this excited state is listed 
in the second column with leading order threoretical uncertainty calculated as described in 

Sec. imn 

Next we looked at a further geometrical property of 2n halo nuclei, specifically the two 
neutron opening angle. As defined in Fig. [6](b), it is straight forward to calculate 6nn using 
the pair of mean square radii found from the three-body wave function \l/c, a/ {r^n) ^^'i 



tan 



Or. 



2 i'^nn) 



(28) 



Our results for known halo nuclei are shown in Table [TTl using the central values of the 
experimental two-body and three-body energies as inputs (compare with Table [T]). We show 
the results for the competing values of the n-c virtual energy for the case of ^^Li. For 
the case of we also show the result with Enc = 0.0 keV, along with its corresponding 
Efimov excited state. The openin g an gle has been calculated from experimental data in two 
recent works by Bertulani et al. 4l[, and Hagino et al. [421], and their results are shown 
in the last two columns, respectively. The study by Bertulani et al. uses the experimental 
values of ■>/ found by Marques et al. j36| seen in the last column of Table [H along 
with two different determinations of a/ ('"c-nn) (see 41] and references within): using laser 
spectroscopy data (results displayed in first row for ^^Li and ^^Be), and also using the B{E1) 
strength (second row for ^^Li). However the result for ^^Be also uses a theoretical calculation 
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IxTi 1 PI 1 1: 


//^2 \ 




ft 

"nn 


/9 [41 1 


"nn I^^J 


11t • 
Li 


8.7±U.7 

^; o_i_i Q 


5.5±(J.4 

n_i_i Q 


77 


roO+lO 

58_M 
66_i8 


r/; 90 + 17.8 

^0-^_21.3 
(;r 90+11. 4 

""J-^-13.0 

[65.29] 




4.1±0.5 


2.8±0.4 


790+16 
'^-13 


64!+^ 






Q n_|_r) a 


1 q_|_n 4 


770+23 
' ' -22 






18c 


2.6±0.7 


1.6±0.4 


78!+f 






20 c 
20c* 


2.8±0.3 
2.7±0.2 
130±10 


2.0±0.2 
l.SibO.l 

69±5 


74!_y 
87!+« 







TABLE II: Various two neutron opening angles of different 2n halo nuclei calculated from the 
results for the mean square radii shown. Compared with results of 41[ and 4^ shown in the last 
two columns, respectively. 



for one of the inputs, rather than being a pure experimental result. On the other hand, 
the study by Hagino et al. uses the experimental values of B{E1) to calculate a/ {rc-nn): 
along with two different determinations of \/Jr^ (see 4^ and references within): using 
experimental values of the nuclear matter radii (first row for ^^Li), and using the results of 
Marques et al. j36| (second row for ^^Li). Also, results found usin g a three-body model with 
density dependant two-body contact interactions were found in [42], and are displayed in 
the third row for ^^Li. As would be expected, our results agree very well with the results 
from the three-body theoretical model used in [45]. Our result for ^^Li using our preferred 
choice of Enc = —800 keV, also agrees very well with the experimental results obtained using 
the value for \/ {^nn)exp Marques et al., which would be expected as our value 
also agrees with this experimental value. Overall there is a good agreement between our 
calculated results and the results of j4l| and 4^ , as all values lie within one error bar of 



each other. However, the size of these error bars suggest that further study should be done 
to improve both the experimental and theoretical results. 



V. CONCLUSION 



In this paper, we have investigated universal aspects of three-body halo nuclei within 
an effective quantum mechanics approach. Assuming that the halo nuclei have resonant 
S-wave interactions between the neutron and the core, the effective potential at leading 
order reduces to a separable S-wave potential. The corrections at next-to- and next-to- 
next-to-leading order in the expansion in Miow I Mugh are determined by the S-wave effective 



ranges [43|, |4^. Corrections from P-wave interactions appear at even higher orders [4:]- An 
important improvement compared to previous calculations is the inclusion of error bands 
based on the omitted higher order terms in the effective theory. 

We have calculated the parametric region within which at least one excited Efimov state 
will occur for different values of the core mass A. The boundary of this region is given by a 
curve in the plane described by the root of the ratio of the two-body bound(virtual) state 
energies to the ground state energy jl9| . We have calculated the boundary of this region for 
various values of the core mass A and provided error bands for the boundary curves. From 
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the current experimental data, we conclude that none of the known halo nuclei is likely to 
have an excited Efimov state. One possible exception is ^''C which could have one excited 
state with a binding energy of less than 7 keV. 

Next, we have studied the structure of known 2n halo nuclei, calculating the one- and 
two-body matter density form factors. From these form factors we were able to extract the 
mean square distances between the two particles in the chosen two-body subsystem, as well 
as the mean square distance of the spectator particle from the center of mass. We found 
that our results for the n-n mean square radius agree well with the experimental data for 
the Borromean halo nuclei ^^Li and ^^Be [36]. We have explicitly not studied the case of 
^He, which is dominated by a P-wave resonance in the n-c interaction ( "^He" ) and requires 
a different counting scheme. While various schemes to treat such P-wave resonances in 



Effective Field Theory have been developed [12, ll3[ , their application to three-body systems 
remains to be worked out. To the expected accuracy, our effective theory gives a good 
description of the studied halo nuclei. Using our results for the mean square distances, we 
have also calculated the two neutron opening angle, and found a good general agreement 



with the recent results of [41i] and |42 



Throughout, this work, we have estimated the theoretical error of the leading order effec- 
tive potential, Eq. ([2]). This uncertainty was quantified in our results through error bands. 
A future study could systematically improve the theoretical error through the inclusion of 
a momentum dependent next-to-leading order term in the effective potential which can be 
matched to the effective range of the interactions. 

Another interesting application of this effective theory will be the study of Coulomb 
excitation data from existing and future facilities with exotic beams (such as FAIR and 
FRIB) . In these experiments a nuclear beam scatters off the Coulomb field of a heavy nucleus. 
Such processes can populate excited states of the projectile which subsequently decay, leading 
to its "Coulomb dissociation" [45]. Effective theories offer a systematic framework for a 
full quantum- mechanical treatment of these reactions. In summary, with new improved 
experimental data for these weakly bound nuclei, much more knowledge can be obtained 
about the structure of these interesting systems as well as discovering whether they show 
universal behavior and excited Efimov states. 
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APPENDIX A: RECONSTRUCTION OF THE WAVE FUNCTION 

In this section we give the expressions for the S-wave part of the full wave functions 
^n{Pi<l) and '^dPiO) where the index i = n,c labels the chosen spectator particle. In 
the wave functions, the p Jacobi momentum describes the relative momentum between the 
two particles in the chosen two-body subsystem, while q describes the momentum of the 
spectator particle relative to the center of mass of the two-body subsystem. 
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For the neutron as the spectator particle, we find: 



1 

1 

Bs)g{nnc)tcKc;B3)F,{n'^,), (Al) 

where the regulator function g and the T-matrices t„ and tc are given in Eqs. ([31 [TT |[T3!) . The 
expressions for the propagators Gq and Gq can be found in Eqs. IQ and ( fTOl) . One can show 
through a simple substitution that Gg (7f„„, 7f^„; Sa) = GQ{p,q; B^), and G'o('7r„c, tt^cj -^s) — 
Gq{p, q; S3), and the propagators simplify. Finally, the shifted momenta are: 



/ N / 1 o A\A + 2)'^ , A(A + 2) , 



<n = q) = \Ip^+ {A + iy ~ aTi ^^^^ 



7r„e = vr„,(p, g) = ^ / - p2 + ^^^^ + ____pq^^ (A4) 



^2 ^ 

<c = KM q) = \ jp'^+ -^.^2 ~ ~aTi^^^^' ^^^^ 

If the core is the spectator particle, we find: 

= j dx GQ{n^n,Kn,B3)g{ncn%{fr'^^-,B3)Fr,{Tr'^J 

q; Bs)g{p)t,iq; B,)FM)- (A6) 

Again, one can show through a simple substitution that Golden, T^'cn'^ B3) = Gq{p, q; -B3). For 
this case, the shifted momenta are given by: 



/ X / A^ , (A + 2y , A(A + 2) 
VTen = vr.„(p, g) = , / ^^^^ p2 + ^^^^^ g2 + ___ 2pg^, (A7) 



TTcn = KniP^ = \l P^ + \q^ ' Pq^- (A8) 
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